Spectra of Linearized Operators for NLS Solitary Waves 
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Abstract. Nonlinear Schrodinger (NLS) equations with focusing power nonlinearities 
have sohtary wave solutions. The spectra of the hnearized operators around these sohtary 
waves are intimately connected to stabihty properties of the sohtary waves, and to the long- 
time dynamics of solutions of (NLS). We study these spectra in detail, both analytically 
and numerically. 
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1 Introduction 

Consider the nonlinear Schrodinger equation (NLS) with focusing power nonlinearity, 

iStV' = -AV' - iV'r" V, (1-1) 

where tp{t, x) : M x — > C and 1 < p < oo. Such equations arise in many physical settings, 
including nonlinear optics, water waves, and quantum physics. Mathematically, nonlinear 
Schrodinger equations with various nonlinearities are studied as basic models of nonlinear 
dispersive phenomena. In this paper, we stick to the case of a pure power nonlinearity for 
the sake of simplicity. 

For a certain range of the power p (see below), the NLS (jl.ip has special solutions, of 
the form tp{t,x) = Q{x)e^^. These are called solitary waves. The aim of this paper is to 
study the spectra of the linearized operators which arise when (jl.ip is linearized around 
solitary waves. The main motivation for this study is that properties of these spectra are 
intimately related to the problem of the stability (orbital and asymptotic) of these solitary 
waves, and to the long-time dynamics of solutions of NLS. 

Let us begin by recalling some well-known facts about (jl.ip . Standard references include 
[H |33l 131] . Many basic results on the linearized operators we study here were proved by 
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Weinstein |381I39| . The Cauchy (initial value) problem for equation (jl.ip is locally (in time) 
well-posed in H^{M."') if 1 < p < Pmaxi where 

Pmax := 1 + 4/(n - 2) if n > 3; Pmax := oo if n = 1, 2. 

Moreover, if 1 < p < pc where 

Pc := 1 + 4/n, 

the problem is globally well-posed. For p > p^, there exist solutions whose ff^-norms go to 
oo {blow up) in finite time. In this paper, the cases p < pc, p = Pc and p > pc are called 
sub-critical, critical, and super- critical, respectively. 

The set of all solutions of (jl.ip is invariant under the symmetries of translation, rotation, 
phase, Galilean transform and scaling: if 'ip{t,x) is a solution, then so is 

ifit, x) := A2/(p-i)V' {XH, XRx - XHv - xq) exp |i 

for any constant xq,v G M", A > 0, 70 G M and R £ 0{n). When p = pc, there is an 
additional symmetry called the "pseudo-conformal transform" (see |34^ p.35]). 
We are interested here in solutions of (jl.ip of the form 

ij{t,x) =Q{x)e'' (1.2) 

where Q{x) must therefore satisfy the nonlinear elliptic equation 

-AQ-\Q\P-^Q = -Q. (1.3) 

Any such solution generates a family of solutions by the above-mentioned symmetries, called 
solitary waves. Solitary waves are special examples of nonlinear bound states, which, roughly 
speaking, are solutions that are spatially localized for all time. More precisely, one could 
define nonlinear bound states to be solutions '4>{t,x) which are non- dispersive in the sense 
that 

sup inf |||x|^/^(t,X - Xo)|| r2mn^ < 00. 

Testing (|1.3p with Q and x.VQ and taking real parts, one arrives at the Pohozaev 
identity ([27\) 

where 

n{p — 1) n + 2— (n — 2)p 

4 4 

The coefficients a and b must be positive, and hence a necessary condition for existence of 
non-trivial solutions is p S (l,Pmax)- 
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For p £ (l,Pmax)) and for all space dimensions, there exists at least one non-trivial 
radial solution Q{x) = Q{\x\) of (|1.3p (existence goes back to H?])- This solution, called 
a nonlinear ground state, is smooth, decreases monotonically as a function of \x\, decays 
exponentially at infinity, and can be taken to be positive: Q(x) > 0. It is the unique 
positive solution. (See [Mj for references for the various existence and uniqueness results 
for various nonlinearities.) The ground state can be obtained as the minimizer of several 
different variational problems. One such result we shall briefly use later is that, for all n > 1 
and p G (l,Pmax); the ground state minimizes the Gagliardo-Nirenberg quotient 

J[u] := ^ jr-^ ^ 1.5 

among nonzero i?i(M") radial functions (Weinstein [55]^. 

For n = 1, the ground state is the unique if^(M)-solution of (jl.Sh up to translation and 
phase [U p. 259, Theorem 8.1.6]. For n > 2, this is not the case: there are countably infinitely 
many radial solutions (still real- valued) , denoted in this paper by Qo,k,pix), /c = 0, 1, 2, 3, . . ., 
each with exactly k positive zeros as a function of |x| (Strauss [32]; see also [2]). In this 
notation, Qo,o,p is the ground state. 

There are also non-radial (and complex-valued) solutions, for example those suggested 
by P. L. Lions [20] with non-zero angular momenta, 

n = 2, Q = (j){r) e*™'^, in polar coordinates r, 9; 

n = 3, Q = </>(r, X3) e*"*^, in cylindrical coordinates r,6,xs, 

and similarly defined for n > 4. When n = 2, some of these solutions are denoted here by 
Qm,k,p with p S (l,Pmax) and k = 0,1,2, .. . denoting their numbers of positive zeros. See 
Section m for more details. 

We will refer to all the solitary waves generated by (5o,o,p as nonlinear ground states, 
and all others as nonlinear excited states. We are not aware of a complete characterization 
of all solutions of ()1.3p . or of (jl.ip . For example, the uniqueness of Qrn,k,p with m, /c ^ 1 
is apparently open. Also, we do not know if there are "breather" solutions, analogous to 
those of the generalized KdV equations. In this paper we will mainly study radial solutions 
(and in particular the ground state), but we will also briefly consider non-radial solutions 
numerically in Section HI 

To study the stability of a solitary wave solution (jl.2p . one considers solutions of (NLS) 
of the form 

tlj{t,x) = [Q{x) + h{t,x)]e'\ (1.6) 

For simplicity, let Q = Qo,o,p be the ground state for the remainder of this introduction (see 
Section |4] for the general case). The perturbation h{t,x) satisfies an equation 

dth = Ch + (nonlinear terms) (1-7) 
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where C is the linearized operator around Q: 



Ch 



i {(-A + 1 - QP-')h - ^ QP-\h + /i)} 



(1.8) 



It is convenient to write £ as a matrix operator acting on 
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(1.9) 



where 



L+ = -A + l-pQ'' 



iP-i 



L_ = -A + 1 - 



iP-i 



(1.10) 



Clearly the operators L_ and play a central role in the stability theory. They are 
self-adjoint Schrodinger operators with continuous spectrum [1, oo), and with finitely many 
eigenvalues below 1. In fact, when Q is the ground state, it is easy to see that L_ is a 
nonnegative operator, while L_|_ has exactly one negative eigenvalue (these facts follow from 
Lemma 12.21 below) . 

Because of its connection to the stability problem, the object of interest to us in this 
paper is the spectrum of the non-self-adjoint operator £. The simplest properties of this 
spectrum are 

1. for all p G (l,Pmax)> is an eigenvalue of £. 

2. the set Sc := {ir : r E M, |r| > 1} is the continuous spectrum of £. 

(See the next section for the first statement. The second is easily checked.) 

It is well-known that the exponent p = pc is critical for stability of the ground state 
solitary wave (as well as for blow-up of solutions). For p < pc the ground state is orbitally 
stable, while for p > Pc it is unstable (see [40^ I12j). These facts have immediate spectral 
counterparts: for p E (l,Pc]) all eigenvalues of £ are purely imaginary, while for p £ 
(PcPmax), £ has at least one eigenvalue with positive real part. 

The goal of this paper is to get a more detailed understanding of the spectrum of £, 
using both analytical and numerical techniques. See [11^ [6l [71 [8] for related work. This finer 
information is essential for understanding the long-time dynamics of solutions of (NLS), for 
example: (i) to prove asymptotic (rather than simply orbital) stability, one often assumes 
either the linearized operator £ has no nonzero eigenvalue, or its nonzero eigenvalues are 
ibri with 0.5 < r < 1. These assumptions need to be verified; (ii) to determine the rate of 
relaxation to stable solitary waves when there is a unique pair of nonzero eigenvalues itri 
with < r < 1. Heuristic arguments suggest that [1/r], the smallest integer no longer that 
1/r, may decide the rate; (iii) to construct stable manifolds of unstable solitary waves, one 
needs to know if there are eigenvalues which are not purely imaginary, and to find their 
locations. These are highly active areas of current research, see e.g. \13\ \T7\ [3T] and the 
references therein. 

Interesting questions with direct relevance to these stability-type problems include: 
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(i) Can one determine (or estimate) the number and locations of the eigenvalues of £ 
lying on the segment between and i? 

(ii) Can ±1, the thresholds of the continuous spectrum Sc, be eigenvalues or resonances? 

(iii) Can eigenvalues be embedded inside the continuous spectrum? 

(iv) Can the linearized operator have eigenvalues with non-zero real and imaginary parts 
(this is already known not to happen for the ground state - see the next section - and 
so we pose this question with excited states in mind). 

(v) Are there bifurcations, as p varies, of pairs of purely imaginary eigenvalues into pairs 
of eigenvalues with non-zero real part (a stability /instability transition)? 

The detailed discussion of the numerical methods is postponed to the Appendix. Roughly 
speaking, we first compute the nonlinear ground state by iteration and renormalization, and 
then compute the spectra of various suitably discretized linear operators. 

Let us now summarize the main results and observations of this paper. 

1. Numerics for spectra. When Q is the ground state, we compute numerically the 
spectra of £, L_|_ and L_ as functions of p, see Figures [THSl In these figures, the 
horizontal axis is the logarithm of p — 1. The vertical axis is: Solid lines are purely 
imaginary eigenvalues of C (without i) for S {l,pc); dashed lines are real eigenvalues 
of C; dotted line are eigenvalues of L^; dashdot lines are eigenvalues of -L_. We have 
ignored imaginary eigenvalues of the discretized operators with modulus greater than 
one, which correspond to the continuous spectra of the original operators. Figured] is 
the one-dimensional case. Figures [2] and [3] are the spectra of these operators restricted 
to radial functions, for space dimensions n = 2 and 3. Figures U] and [5] are for n = 2 and 
are the spectra restricted to functions of the form 0(r)e*^ and (j){r)e^'^^ , respectively. 
These pictures shed some light on questions (i), (iv), and (v) above, and to a certain 
extent on question (ii). 

Figures [T0HT5] are concerned with the spectra of excited states, see discussion below. 

2. One-dimensional phenomena. The case n = 1 is the easiest case to handle ana- 
lytically. In Section [3l we undertake a detailed study of the one-dimensional problem, 
giving rigorous proofs of a number of phenomena observed in Figure [H One simple 
such phenomenon is the (actually classical) fact that the eigenvalues of and L_ 
exactly coincide, with the exception of the first, negative, eigenvalue of L-|_ (note that 
this appears to be a strictly one-dimensional phenomenon: the eigenvalues of L_|_ and 
L_ are different for n > 2, as Figures [JHS] indicate) . In fact, we are able to prove 
sufficiently precise upper and lower bounds on the eigenvalues of £ (lying outside 
the continuous spectrum) to determine their number, and estimate their positions, 
as functions of p (see Theorem 13. Sh . We use two basic techniques: an embedding of 

and -L_ into a hierarchy of related operators, and a novel variational problem for 
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p 

1 .04 1.1 1 .32 1 .82 2 3 5 




-1.39 -1 -0.5 -0.086 O 0.301 0.5 0.65 



logiob - 1) 

Figure 1: Spectra of C, and L_ for n = 1 with logarithmic axis for the values of p — 1. 
(solid line: purely imaginary eigenvalues of £; dashed line: real eigenvalues of C\ dotted 
line: eigenvalues of L_|_; dashdot line: eigenvalues of L_) 



P 

1.01 1.03 1.1 1.32 2 2.27 2.379 3 




-2 -1.5 -1 -0.5 0.34 

logio(p- 1) 

Figure 2: Spectra of C, L_|_ and L_ restricted to radial functions in the two-dimensional 
space, with logarithmic axis for the values oip—1. (solid line: purely imaginary eigenvalues 
of C\ dashed line: real eigenvalues of C\ dotted line: eigenvalues of L_)_; dashdot line: 
eigenvalues of L_) 
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p 

1 .01 1 .03 1 .1 1 .32 3.24 2 3.19 




-2 -1 .5 -1 -0.5 -0.3497 0.34 



logio(p - 1) 

Figure 5: Spectra of £, and L_ restricted to functions of the form (/)(r)e*^^ in the 
two-dimensional space. 

the eigenvalues, in terms of a 4-th order self-adjoint differential operator (see The- 
orem I3.6p . In this way, we get a fairly complete answer to question (i) above for 
n = \. 

3. Variational characterization of eigenvalues. We present self-adjoint variational 
formulations of the eigenvalue problem for C in any dimension (see Summary 12. Sp , 
including the novel n = 1 formulation mentioned above. In principle, these provide 
a means of counting/estimating the eigenvalues of C (and hence addressing question 
(i) above in higher dimensions), though we only obtain detailed such information for 
ra = 1. 

4. Bifurcation at p = pc. In each of Figures [THHl a pair of purely imaginary eigenvalues 
for p < Pc appears to collide at at p = pc, and become a pair of real eigenvalues for 
p > Pc- This is exactly the stability/instability transition for the ground state. We 
rigorously verify this picture, determining analytically the spectrum near for p near 
Pc, and making concrete a bifurcation picture suggested by M. I. Weinstein (personal 
communication): see Theorem 12.61 This gives a partial answer to question (v) above. 
It is worth pointing out that for n = 1, the imaginary part of the (purely imaginary) 
eigenvalue bifurcating for p < pc is always larger than the third eigenvalue of (the 
first is negative and the second is zero) - this is proved analytically in Theorem 13.81 
For n > 2, however, they intersect at p ~ 2.379 for 2D and p ~ 2.046 for 3D, (see 
Figures HHSD. 
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5. Interlacing property. A numerical observation: in all the figures, the adjacent 
eigenvalues of L seem each to bound an eigenvalue of and one of -L_ (at least for 
p small enough). We are able to establish this "interlacing" property analytically in 
dimension one (see Theorem I3.8p . 

6. Threshold resonance. An interesting fact observed numerically (Figured]) is that, 
in the ID case, as p ^ 3, one eigenvalue curve converges to iti, the threshold of 
the continuous spectrum. One might suspect that, at p = 3, ±z corresponds to a 
resonance or embedded eigenvalue. It is indeed a resonance: we find an explicit 
non spatially-decaying "eigenfunction" , and show numerically in Section 13.71 that the 
corresponding eigenfunctions converges, as p ^ 3, to this function. This observation 
addresses question (ii) above for n = 1. 

7. Excited states. In Section [J] we consider the spectra of linearized operators around 
excited states with non-zero angular momenta. We observe that, in addition to the 
bifurcation mentioned above at p = pc, there are complex eigenvalues which are nei- 
ther real nor purely imaginary (addressing question (iv) above; see Figures [T0Hl5]l . 
symmetric with respect to both real and purely imaginary axes. These complex eigen- 
values also come from bifurcation: as p decreases, a quadruple of complex eigenvalues 
will collide into the imaginary axis away from 0, and then split to 4 purely imaginary 
eigenvalues. It seems that all eigenvalues lie on the imaginary axis for p E (l,p*) 
for some p=K close to 1. In other words, numerically these excited states are spectrally 
stable for p close to 1. It is possible that the numerical error increases enormously as 
p ^ 1-1- due to the artificial boundary condition, since the spectrum is approaching to 
the continuous one for p = 1. This has to be verified analytically in the future. Even 
if they are indeed spectrally stable, it is not clear if they are nonlinearly stable. 

It is worth mentioning some important questions we cannot answer: 

1. We are so far unable to give precise rigorous estimates on the number and positions 
of the eigenvalues of £ for n > 2 (question (i) above). 

2. We cannot exclude the possible existence of embedded eigenvalues (question (iii) 
above). 

3. We do not know a nice variational formulation for eigenvalues of C when Q is an 
excited state (this problem is also linked to question (i) above). 

4. We do not have a complete characterization of solitary waves, or more generally of 
nonlinear bound states. 

We end this introduction by describing some related numerical work. Buslaev-Grikurov 
O study the linearized operators for solitary waves of the following ID NLS with p < q, 
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They draw the bifurcation picture for eigenvalues near zero when the parameter a > is 
near a critical value, with the frequency of the solitary wave fixed. This picture is similar 
to Weinstein's picture which we study in Section [2.31 

Demanet and Schlag [8] consider the same linearization as us and study the super-critical 
case n = 3 and p < 3 near 3. In this case, it is numerically shown that both and L_ 
have no eigenvalues in (0, 1] and no resonance at 1, a condition which implies (see [3l]) that 
C has no purely imaginary eigenvalues in [— i,0) U (0, and no resonance at ±1. 

We outline the rest of the paper: in Section [2] we consider general results for all dimen- 
sions for ground states. In Section [3] we consider one dimensional theory. In Section H] we 
discuss the spectra for excited states with angular momenta. In the Appendix we discuss 
the numerical methods. 

Notation: For an operator A, N{A) = |(/) G A(j) = O} denotes the nullspace of A. 
Ng{A) = ^'^iN{A^) denotes the generalized nullspace of A. The L^-inner product in 
is {f,g) = f^n fgdx. 

2 Revisiting the general theory for ground states 

In this section we review mostly well-known results which are valid for all dimensions, for 
the ground state Q{x) = Qo,o,p{x), and give new proofs of some statements. 

We begin by recalling some well-known results for the linearized operator £ defined by 
(jl.Sp . As is well known for linearized Hamiltonian system (and can be checked directly) , if A 
is an eigenvalue, then so are —A and itA. Hence if A 7^ is real or purely imaginary, it comes 
in a pair. If it is complex with nonzero real and imaginary parts, it comes in a quadruple. It 
follows from nonlinear stability and instability results [401 12] that all eigenvalues are purely 
imaginary if p £ {l,pc), and that there is at least one eigenvalue with positive real part 
when p £ (PcPraax)- It is also known (see e.g. [7]) that the set of isolated and embedded 
eigenvalues is finite, and the dimensions of the corresponding generalized eigenspaces are 
finite. 

2.1 L+, L , and the generalized nullspace of C 

Here we recall the makeup of the generalized nullspace Ng[C) of C. Easy computations give 
L+Qi = -2Q, L_Q = 0, where Qi := (^ + X • V)Q, (2.1) 

and 

L^xQ = -2VQ, L+VQ = 0. (2.2) 

In the critical case p = Pc, we also have 

L_(|a;pQ) = -4Qi, L+p=\x\'^Q (2.3) 

for some radial function p{x) (for which we do not know an explicit formula in terms of Q). 
Denote 
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For 1 < p < Pmax J the generalized nullspace of £ is given by (see [39 
Ng[C) = span 





xQ 





\x?Q 



VQ 




Qi 




6P 

Pc 



(2.5) 



In particular 

dim Ng{C) = 2n + 2 + 25l^. 

The fact that the vectors on the r.h.s of (j2.5p lie in Ng{C) follows immediately from the 
computations ()2.ip - p.3p . That these vectors span Ng{C) is established in [39], Theorems 
B.2 and B.3, which rely on the non-degeneracy of the kernel of L^: 



Lemma 2.1 For all n> 1 and p G {l,pma.x), 

N{L+) =span{VQ} 

This lemma is proved in [39] for certain n and p (n = 1 and l<p<oo, orn = 3 and 
1 < p < 3), and is completely proved later by a general result of p^. We present here a 
direct proof of this lemma, without referring to [18] , relying only on oscillation properties of 
Sturm-Liouville ODE eigenvalue problems. A similar argument (which in the present case, 
however, applies only for p < 3) appears in [H], Appendix C. For completeness, we also 
include some arguments of [39]. 

A new proof. We begin with the cases n>2. Since the potential in Lj^ is radial, any 
solution of Lj^v = can be decomposed as v = l^fc>o SjeEfe where r = 

X = ^ is the spherical variable, and Y/jj denote spherical harmonics: — Agn-iYfej = Afclfej 
(a secondary multi-index j, appropriate to the dimension, runs over a finite set for each 
k). Then L+f = can be written as A^f/j j = 0, where, for = 0, 1, 2, 3, . . ., 

Ao = -dl-'^^—^dr + l-pQP-^{r), Ak = Aq + Xkr-"^ , Xk = k{k + n-2). 
r 

Case 1: k = \. Note VQ = Q'{r)x. Since AiQ' = and Q'{r) < (monotonicity of the 
ground state) for r £ (0, oo), Q'{r) is the unique ground state of Ai (up to a factor), and 
so Ai > 0, ^i||Q,}± > 0. 

Case 2: k > 2. Since Ak = Ai + {Xk — Ai)r~^ and > Ai, we have Ak > 0, and hence 
AkVk = has no nonzero L^-solution. 

Case 3: k = 0. Note that the first eigenvalue of is negative because {Q, AqQ) = 
(Q) ~iP ~ ^)Q^) < 0. The second eigenvalue is non-negative due to ()2.7p and the minimax 
principle. Hence, if there is a nonzero solution of ^40^0 = 0, then is the second eigenvalue. 
By Sturm-Liouville theory, vo{r) can be taken to have only one positive zero, which we 
denote by tq > 0. By (gl]), AqQ = -{p - 1)Qp and A^Qi = -2Q. Hence (Qp,wo) = = 
{Q,vq). Let a = (Q(ro))P"^ Since Q'(r) < for r > 0, the function QP-aQ = Q{QP-^-a) 
is positive for r < tq and negative for r > tq. Thus vq{QP — aQ) does not change sign, 
contradicting (vq, Qp — aQ) = 0. Combining all these cases gives Lemma |2. II for n > 2. 
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Finally, consider n = 1. Suppose L+v = 0. Since preserves oddness and evenness, 
we may assume v is either odd or even. If it is odd, it vanishes at the origin, and so by 
linear ODE uniqueness, f is a multiple of Q'. So suppose v is even. As in Case 3 above, 
since has precisely one negative eigenvalue, and has Q' in its kernel, v{x) can be taken 
to have two zeros, at x = ±xq, xq / 0. The argument of Case 3 above then applies on 
[0, oo) to yield a contradiction. □ 

We complete this section by summarizing some positivity estimates for the operators 
and L_. These estimates are closely related to the stability/instability of the ground 
state. 



Lemma 2.2 



max J 

(Q,L+Q)<0, L+\^Qpy±>0 (l<p<Pmax), 

L+|,^,x>0 {l<p<pc) 



L. 



-^{Q,xQ} 
{Q,xQ,\x\2Q} 



> 0, L_ 
> 0, 



{Qi}->0 (l<P<Pc) 



(2.6) 
(2.7) 
(2.8) 
(2.9) 
(2.10) 



Proof. Most estimates here are proved in [39] except the second part of (j2.7p when p > pc- 
It can be proved for p £ (l,Pmax) by modifying the proof of [391 Prop. 2.7] for ()2.8I) as 
follows. (It is probably also well-known but we do not know a reference.) 

Recall the ground state Q is obtained by the minimization problem (II. 5p . If a minimizer 
Q{x) is rescaled so that 



I\VQ\' _IQ' _JQ 



p+i 



2a 



2b 



P+l 



constant A; > 0, 



i.e., (II. 4p is satisfied, then Q{x) satisfies ()1.3p . The minimization inequality ^|^_ 
£7]] > for all real functions rj, is equivalent to 

kin, L+v) >\j '?Ag)2 + Qrif - ( j Q^rif. 



(2.11) 



Thus (?7, L+ri) > if _L Q^. Note that, if ?7 _L Q, by ([L3]) the right side of (pTTT) is 
positive if a < 1, i.e. p < Pc- In this way, we recover ()2.8p . □ 



2.2 Variational formulations of the eigenvalue problem for C 

In this subsection we summarize various variational formulations for eigenvalues of C. The 
generalized nullspace is given by (j2.5p . Suppose A 7^ is a (complex) eigenvalue of £ with 
corresponding eigenfunction [^] G , 



(2.12) 






LJ 




u 


= A 


u 















w 


w 
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The functions u and w satisfy 

L+n = —Xw, L^w = An. (2-13) 

Therefore 

L^L^u = fiu, ij, = -\^. (2.14) 

Since (/xu, Q) = {L^L^u, Q) = {L^u, L^Q) = and fi ^ 0, we have u ± Q. 

Denote by 11 the L^-orthogonal projection onto Q^. We can write L^u = HL^u + aQ. 
Eq. (|2.14p implies L^HL^u = jju and hence, using u -L Q and ()2.6p . IlL^u = Ll^/xu. Thus 

{u,Q) = 0, L+u = fiLZ^u + aQ. (2.15) 

Since (|2.14p is also implied by Eq. (j2.15p . these two equations are equivalent. 

If Q{x) is a general solution of (jl.Sp . /x = — may not be real. However, it must be 
real for the nonlinear ground state Q = Qco.p- This fact is already known (see [29J). We 
will give a different proof. 

Lemma. For Q = Qo,o,p^ every eigenvalue /i of (j2.14p is real. 

A new proof. Multiply (j2.13p by u and w respectively and integrate. Then we get 



{u, Lj^u) = —X{u,w), {w, L-w) = X{w,u) = X{u,w). (2-16) 

Taking the product, we get 

{u, Lj^u){w, L^w) = — A^|(n,zi;)p = fi\{u,w)\'^ . 

If /i 7^ 0, If is not a multiple of Q, and so by (|2.6p . {w,L^w) > 0. Hence {u,w) / by 
dZin]). Thus 

{u,L+u){w,L_w) 

M - r? ^ ^• 

|(m,w)|^ 

□ 

This argument does not work when Q is an excited state, since (n, w) may be zero (see 
e.g. [371 Eq.(2.63)]). The fact n £U. implies that eigenvalues A of £ are either real or purely 
imaginary. Thus C has no complex eigenvalues with nonzero real and imaginary parts. This 
is not the case for excited states (see Section HJ also ^37j ) . 

The proof of reality of /i in [29] uses the following formulation. For the nonlinear ground 
state Q, L_ is nonnegative and the operator L^J'^ is defined on and invertible on Q-^. 
A nonzero /x G C is an eigenvalue of (j2.14p if and only if it is also an eigenvalue of the 
following problem: 

L^J^L+L'J^g = i^g, (2.17) 

with g = L_^^'^u. The operator L^/^L+Li^^ already appeared in [36]. Since it can be 
realized self-adjoint operator, must be real. 
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1/2 1/2 

Furthermore, the eigenvalues of L_ L^L_ can be counted using the minimax princi- 
ple. Note that Q is an eigenfunction with eigenvalue 0. For easy comparison with other 
formulations, we formulate the principle on Q^. Let 

/x,:= inf j = l,2,3,... 2.18 

<;±Q,gfc,fc=l,...j-l 

with a suitably normalized minimizer denoted by gj (if it exists - the definition terminates 
once a minimizer fails to exist). The corresponding definition for (I2.15[) is 

^, := mf ) ' _+ \ , (j = 1,2,3,...) (2.19) 

M±Q,(u, LI^Ufc)=0, fe=l,...,j-l [U,L_ U) 

with a suitably normalized minimizer denoted by Uj (if it exists). In fact, the minimizer Uj 
satisfies 

L-^-Uj = ^jLZ^Uj + ajQ + PiLZ^ui + • • • + Pj-iLZ^Uj-i, (2.20) 

for some Lagrange multipliers . . . Testing ()2.20p with Uk with k < j, we get 
{uk, PkLZ^Uk) = {uk,L-^-Uj) = {L-\.Uk,Uj) = by (I2.20p for Uk and the orthogonality condi- 
tions. Thus = and L+Uj = fijLZ^Uj + ajQ and hence uj satisfies (j2.15p . 

Lemma 2.3 The eigenvalues of (j2.18p and (j2.19p are the same, and 



ifl<p<Pc 
ifp = Pc 

ifPc<P< Pmax 



/il = • • • = /in = 0, > 0. 

Hl = ■■■ = fJ-n+1 = 0, /i„,+2 > 0. 

/"I < 0, /U2 = • • • = /Xn+1 = 0, Hn+2 > 0. 



T/ie 0-eigenspaces are spanL_^^'^{VQ,Sp^Qi} for (j2.18p and span{VQ,6p^Qi} for (j2.19p , 
where 5p^ is defined in ()2.4p . 

Proof. The eigenvalues of (j2.18p and (j2.19p are seen to be the same by taking g = 
L_ w up to a factor. By estimate (|2.8p . /ii > for p G (l,Pc]- For p G (PoPmax); using 
(|1.4p . nQi = Qi — ^(^^Q) ^-ud elementary computations (such as (j2.22p below), one finds 

(nQi,L+nQi) = ""'^P - ^\ p,-p)^ [ Qp+^ 

4 P+l J 

which is negative for p > pc- Thus fii < 0. By estimate ()2.7p . /i2 > for p G (l,Pmax)- 
It is clear that u = gf-Q, j = 1, . . . , n, provides n 0-eigenfunctions. For p = Pc, another 

0-eigenfunction is u = Qi since Qi -L Q (see again (j2.22p below), LI^VQ = —^xQ, and 
{Qi, L-^-Qi) = 0. It remains to show that > for p G (l,Pc) and /i„,+2 > for 

p G [pcjPmax)- If fJ'n+2 = for p G (PcPmax), the argument after ()2.20p shows the existence 
of a function Un+2 7^ satisfying 

= for some a G M, u„+2 -L Q, LZ^ui, LZ^VQ = —\xQ. 

By Lemma |2.H Un+2 + f Qi = c • VQ for some c G M*^. The orthogonality conditions imply 
Un+2 = 0. The cases p G (l,Pc] are proved similarly. □ 
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Remark 2.4 The formulation (j2.19p for /Ui has been used for the stabihty problem, see 
e.g. |34l P-73, (4.1.9)], which can be used to prove that //i < if and only if p G (PciPmax) 
by a different argument. The later fact also follows from [69\ \n\ indirectly. 

We summarize our previous discussion in the following theorem. 

Summary 2.5 Let Q{x) he the unique positive radial ground state solution of (jl.3p . and 
let C, and L_ he as in (|1.8|) and (|1.1U|) . The eigenvalue prohlems (|2.14|) . (|2.15|) . and 
(j2.17p for 7^ are equivalent, and the eigenvalues fi must he real. These eigenvalues can 
he counted hy either (I2.18P or (12.19p . < if and only if p £ (pcPmax)- Furthermore, all 
eigenvalues of £ are purely imaginary except for an additional real pair when p G (pciPmax)* 

The last statement follows from the relation /i = — in ()2.14p . 



2.3 Spectrum near for p near pc 

We now consider eigenvalues of C near when p is near pc- It was suggested by M.I. 
Weinstein that as p approaches pc from below, a pair of purely imaginary eigenvalues will 
collide at the origin, and split into a pair of real eigenvalues for p > pc- In the following 
theorem and corollary we prove this picture rigorously and identify the leading terms of the 
eigenvalues and eigenfunctions. 

Note that Comech-Pelinovsky [6] considers a different problem where the equation is 
fixed and the varying parameter is frequence uj rather than exponent p of the nonlinearity. 
That problem has only U{1) symmetry and no translation, but its situation is similar to 
ours since we consider radial functions only in our proof. It seems one can adapt their 
approach to give an alternative proof. They use an abstract projection (Riesz projection) 
onto the discrete spectrum to reduce the problem to a 4x4 matrix problem (and exploit the 
complex structure), while we are more direct. We thank the referee for pointing out [6] to 
us. 

Theorem 2.6 There are small constants /i* > and > so that for every p G {pc — 
£*,Pc + there is a solution of 

L+L^w = fiw (2.21) 

of the form 

w = Wo + {p - Pcfg, Wo = Q + a{p - pc)\x\'^Q, g-\-Q, 
= 8a{p - Pc) + {p - Pcfv, a = a{p) = ^q^'^q-^ < 0^ 

with ||g'||j;^2, \'r]\, \a\ and l/\a\ uniformly hounded in p. Moreover, for p ^ pc, this is the 
unique solution of ()2.2ip with < |^| < 

Proof. Set e : = p — pc- Computations yield 

iQuQ) = - ^) (Q, Q) = -^(^(^' (2.22) 
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(Qi, Q") = ^—r{L+Q, Qi) = ^—{Q, L+Qi) = -^{Q, Q), (2.23) 

p — 1 p — 1 p — 1 

and 

(Qi, \x\'Q) = (^-^ - ^) (Q, |x|2g) = -(1 + ^^^)iQ, \^\'Q)- (2.24) 

Since by (l2:2TD with /i / 0, 

iQi,w) = fi^^{Qi,L+L^w) = fi^^{L-L+Qi,w) = 0, 

we require the leading term {Qi,wo) = 0, which decides the value of a using (j2.22p and 
(j2.24p . Thus we also need {Qi,g) = 0. That a < (at least for e sufficiently small) follows 
from ([2:23]) and ([221]). 
Using the computations 

L_|xPQ = [L_, IxPlQ = -4x ■ VQ - 2nQ = -4Qi - -^^eQ (2.25) 

p — 1 

and 

L+Q = [L_ -{p- 1)QP-']Q = -{p- 1)QP, 

we find 

2n 

L+L^wq = aeL+l-AQi eQ] = ae[8Q + 2neQP]. 

p — 1 

Thus fj, = 8ae + o(e) and we need to solve 

= [L+L_ - 8ae - e'^r]][wo + e'^g] 

which yields our main equation for g and r]: 

L+L^g = 8a^{\x\'^Q) - 2an{QP) + r]Wo + (8ae + e^r])g. (2.26) 

Recall that on radial functions (we will only work on radial functions here) 

ker[(L+L_)*] = ker[L_L+] = span{Qi}. 

Let P denote the L^-orthogonal projection onto Qi and P := 1 — P. It is necessary that 

P[8a'^{\x\'^Q) - 2an{QP) + r]Wo + (8ae + e'^r])g] = 

for (I2.26P to be solvable. This solvability condition holds since {Qi,g) = (QijWq) = 0, and, 
using the relations (|2:2iD and (|2:23D . (Qi, Sa^dxpQ) - 2an{QP)) = 0. 
Consider the restriction (on radial functions) 

T = L+L_ : [kerL_]-L = — > Ran{P) = Qf. 
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Its inverse T ^ = ""^(1^+) ^ is bounded because ^ : Qi — > Q"*" and ^ : 

Q-*- —>■ are bounded. So our strategy is to solve ()2.26p as 

g = r~^P[8a2(|xpQ) - 2an(QP) + r]WQ + (8ae + e'^ri)g\ (2.27) 

by a contraction mapping argument, with ry chosen so that {Qi,g) = 0. Specifically, we 
define a sequence 5o = 0) % = 0, and 

5fc+i = PT-^P[8aWx\^Q) - 2an{QP) + %w;o + (8ae + e^^)^^], 

m+i = -TTT^^T-^ {Qi,T-^P[8a\\x\^Q) - 2an{QP) + (Sae + e^r]k)gk])- 
{Qi,T ^wo) 

We need to check {Qi,T~^wo) is of order one. Since wq = Q + 0{e) and L+Qi = —2Q, 
we have {L^)~^wq = — ^IIQi + 0(e) where 11 denotes the orthogonal projection onto Q-^. 
Thus, using (f2:25|l and ([2221), 

{Qi,T-^wo) = -^{Qi, (L^r^UQi) + 0{e) 

= i(Qi,n|x|2Q) + 0(£) = l(Qi,|x|2Q)+0(e), 

which is of order one because of (j2.24p . One may then check that := ||5fc+i — 5fc|lL2 + 

£1/2 1 

Vk+i - Vk\ satisfies A^'^+i < Ce^/^iVfc, and hence {gk,rik) is indeed a Cauchy sequence. 
Finally, the uniqueness follows from the invariance of the total dimension of generalized 
eigenspaces near under perturbations. □ 

Remark 2.7 To understand heuristically the leading terms in w and fi, consider the fol- 
lowing analogy. Let = || ^ , which corresponds to L^L^. One has [J] = [[]], 

[5] = [g] and [\] = £[\]. The vectors [q], [\] and [\] correspond to Q, \x\'^Q and w, 
respectively. 

The theorem yields an eigenvalue ^ with the same sign as pc — p- Since the eigenvalues 
of C are given by A = ziz^—fi, we have the following corollary. 

Corollary 2.8 With notations as in Theorem \2.6l C has a pair of eigenvalues A = ib-^— ^ = 
±iy8|a|(p — Pc) — (p — Pc^rj with corresponding eigenvectors [^] solving (j2.12p and 

u = X-^L^w = ^^/2\a\{p-pc)Qi + 0{{p - pc)^/^). 

When p £ {pc — £*, Pc) (stable case), A and u are purely imaginary. 
When p £ (pc, Pc + £*) (unstable case), A and u are real. 

In deriving the leading term of u we have used (I2.25p . We solved for w before u simply 
because w is larger than u. 
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3 One dimensional theory 

In this section we focus on the one dimensional theory. For n = 1, the ground state Q{x) 
has an exphcit formula for all p G (1, co), 

Q(x) = CpCOsh-^(x//3), cp:={^)^, (3:=^. (3.1) 

2 p — 1 

The function Q{x) satisfies (jl.3p and is the unique ff^(M)-solution of (jl.Sp up to translation 
and phase \^ p. 259, Theorem 8.1.6]. 

3.1 Eigenfunctions of and L_ 

We first consider eigenvalues and eigenfunctions of Lj^ and L_. For n = 1, 

L+ = -a,, + 1 - pQP-\ L_ = -a,, + 1 - QP~\ (3.2) 

By (13. ip . these operators are both of the form 

-a^^ + l-Csech2(x//3). 

Such operators have essential spectrum [l,oo), and finitely many eigenvalues below 1. A 
lot of information about such operators is available in the classical book [35j, p. 103: 

• all eigenvalues are simple, and can be computed explicitly, as zeros and poles of an 
explicit meromorphic function; 

• all eigenfunctions can be expressed in terms of the hypergeometric function. 

We begin by presenting another way to derive the eigenvalues, as well as different for- 
mulas for the eigenfunctions. We will not prove right here that this set contains all of the 
eigenvalues/eigenfunctions. This fact is a consequence of the more general Theorem 13.41 
proved later (and see also [35]). 

Define 

A - 1 - A-2 k - ^±i- - "^^P-^) 

^1 2 2 (3 3) 

m + 1 „ ^ 
Pm '■= tor m > 1, pi = Qo. 

m — 1 

The following theorem agrees with the numerical observation Figure [TJ 

Theorem 3.1 For n = 1 and 1 < p < oo, let Q{x) he defined by ()3.ip . and L- be 
defined by ()3.2p . and Xmi^miVm be defined by ()3.3p . Suppose for M G 

VM+i <P<PM- (3.4) 



18 



Then the operator Lj^ has eigenvalues Am, < m < M , with eigenf unctions of the form 

e e 

j=0 j=l 
and the operator L_ has eigenvalues \m, 1 < m < M, with eigenf unctions of the form 

e. e 

In particular, all eigenvalues of L_ are eigenvalues of L+, and always has one more 
eigenvalue (Xq <0) than L_. 

Proof. It can be proved by induction, using 

p_i _ 2/^/r,^ r) - -n+^-nhf'r/R^ rP■-rP■(^ ?_np-i^ 



= cosh~^(x//3), = -Qtanh(x//3), Qi = Q'{1 -Q 

2 p+l 



and 



p+l 
p+l 

Q-kj^ Qk ^ {k-mk+p+i) ^p„i ^ ^2 

p + l 

The coefficients of Q^~^ vanish when /c = 1^ 1^ 2z£^ respectively. It is why the highest 

p+l 3— p 

power of Q is Q 2 in (/?2£, Qx in </'2£-i, Q in V'2£-i, and (Q 2 j^. in V'2f- □ 



3.2 Connection between L+ and L and their factorizations 

In hght of Theorem 13. H it is natural to ask why all eigenvalues of L_ are also eigenvalues 
of L+. Is there a simple connection between their eigenfunctions? In this section we prove 
this is indeed so. 

We first look for an operator U of the form 

U = d^ + R{x), {soU* = -d^ + R{x)), 

such that 

L_U = UL+, { so U* = L+U*). (3.7) 
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It turns out that there is a unique choice of R{x): 

p+lQx {p-l)x 
R{-) = = ^tanh(^^). 

In fact, with this choice of R{x), 

U = (/?o5a;93o ^ ( so f7* = -(/?o ^9^(y9o), (3.8) 

where 'fQ = Q^ is the ground state of L+, and is considered here as a multiphcation 
operator: U f = ipod^i^Q^ f) ■ 

Suppose now ■i/' is an eigenfunction of L_ with eigenvalue A: L-^j: = Xip. By (j3.7p . 

= C/*(L_ - A)^ = (L+ - X)U*ip. 

Thus U*^p is an eigenfunction of with same eigenvalue A (provided U*Tp G L^). Therefore, 
the map 

^ U*ip 

sends an eigenfunction of L_ to an eigenfunction of L+ with same eigenvalue. This map is 
not onto because U* is not invertible. Specifically, the ground state ipo is not in the range. 
In fact, Uifo = ifodxifQ^ipo = 0. If ipo = U*ijj, then {ipo,ipo) = {ipo,U*ip) = {Uipo,ip) = 0, a 
contradiction. We summarize our finding as the following proposition. 

Proposition 3.2 Under the same assumptions and notation as Theorem \3.1l the eigen- 
functions ipm and ipm of L+ and L_ satisfy 

iPm = U*'lpm, (m = l,...,M), 

up to constant factors. Note that U* sends even functions to odd functions and vice versa. 

Proof. We only need to verify that U*iprn G L^- This is the case since U* = —dx + 
2ii tanh(x//3), ipm{x) are sums of powers of Q and Qz, and that tanh(x//5), Qx/Q, and 
Qxx/Qx are bounded. □ 

Analogous to the definition of U , we define 

S ■.= QdxQ-^ = dx-^, {so S* = -Q-^dxQ). (3.9) 

Clearly SQ = 0. Recall that Aq is the first eigenvalue of L-|_ with eigenfunction ipQ. Hence 
— Aq is a nonnegative operator. In fact we have the following factorizations. 

Lemma 3.3 Let U and S be defined by (|3.8p and (|3.9p . respectively. One has 

L+-Xo = U*U, L^-Xo = UU*. (3.10) 

L_=S*S, SS* = -dl + l + ^^^QP-\ (3.11) 

Moreover, SS* = L_ + ^^^Q^-^ > 0. 
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The formula L_ = S*S was known, see e.g. |34^ p. 73, (4.1.8)]. It is an example of 
the Darboux transformations, see e.g. [22]. Factorization of Schrodinger operators into 
first-order operators has been known since the times of Darboux (1840s). 



3.3 Hierarchy of Operators 

In this subsection we generahze Theorem 13.11 and Lemma 13.31 to a family of operators 
containing and -L_. As a reminder, we have 



Q"/Q = l-QP-\ (Q'/Qf = l-^QP-\ 

p+l 

iQ'/QY = Q"/Q - {Q'/Q? = -^Q^'^ 



More explicitly. 



(3.12) 



Let S{a) := Q^d^Q-". We have 

S{a) = dx- aQ'/Q, S{a)* = -d^ - aQ'/Q, 

SiaySia) = -a^ + a2 - a |a + ^1 ^Q^"'- ^^'^^^ 
Define the fohowing hierarchy of operators: 

Sj := S{kj), where recall kj = 1 — (j — 1)^ 



2 (3.14) 

Lj := Sj-iS*_i + Xj-i = S*Sj + Xj, where recall Xj = 1 — kj. 

Then we have 



5*0 — U, Si — S, . . . 

Lo = L+, Li = L^, L2 = SS*,... (3.15) 



Lj = -dl + l- k,.ik,-^QP-\ (3.16) 



Note that j here can be any real number. 

Recall the definition pj := 1 + 2/(j — 1) for j > 1, and set pj = oo for j < 1. Then pj 
is a monotone decreasing function of j, kj > for p < Pj, kj = for p = pj and kj < for 
p > Pj- Let 

Xj {l<P<Pj), 

1 {pj<p<pj-i), (3.17) 
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By the second identity of (j3.14p . and (j3.16p together with the fact kj^ikj < for pj < p < 
Pj-i, we have the lower bound 

Lj > X'j. (3.18) 

In fact, this estimate is sharp: for p G U (j)j_i,oo), the ground state is obvious from 

the second identity of (j3.14p : 

(ijQj = XjQj, {l<p< Pj), 

where we denote 

Qj-=Q'"', Q*:=Q-^'. (3.20) 

For p G \pj,pj-i\, there is no ground state. Thus we have completely determined the ground 
state of Lj for all p > 1. The complete spectrum, together with explicit eigenfunctions, are 
derived using the third identity of (j3.15p as follows. 

Theorem 3.4 For any j G M and p > 1, the point spectrum of Lj consists of simple 
eigenvalues 

speCp{Lj) ={\k \ p<Pk, {j,j + 1, j + 2, . . . }} 

U {Xk \ p > Pk, kG {j - 1, j - 2, j - 3, . . .}}, 



and the eigenfunction for the eigenvalue Xk is given uniquely up to constant multiple by 

(3.22) 



S*j---Sl^,Qk {ke{j,j + l,...}), 
S,-i---Sk+iQl (feG{j-l,j-2,...}), 



each of which is a linear combination of 



Qj,Qj+2,---Qk {k e + 2, ...}), 

Qj+iR, Qj+3R, ■ ■ ■ QkR (k G {i + 1, j + 3, . . . }), 

Q*j.i, Q*j.3^ ...Ql {ke {j - 1, j - 3, . . . }), 

[Q*_^R,Q*_^R, . . . QIR {k G {j - 2,j - 4, . . . }) 



(3.23) 



where R := Q' /Q. 



Proof. The ground states have been determined. The third identity of (j3.15p implies 
that ()3.22p belong to the eigenspace of Lj with eigenvalue Xk- Moreover, each function is 
nonzero because is injective for p < Pk and so is Sk for < p. Since Sj annihilates only 
the ground state Qj for p < pj and 5'*_x annihilates only the ground state Q*j-i for p > Pj, 
all the excited states of Lj for p < pj are mapped injectively by Sj to bound states of Lj+i, 
and for p > pj by Sj_i to those of ij-i. Hence we have (I3.2ip and all the eigenvalues are 
simple because the ground states are so. ()3.23p follows from the fact that Sj and S* act on 
Q" like C{a,j)R, while SjSj-i and ^^.i^* act on Q'' like Ci{a,j) + C2{aJ)QP-^. □ 
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3.4 Mirror conjugate identity 

The following remarkable identity has application to estimating eigenvalues of L (see Sec- 
tion [M]): 

- A,)5* = 5;(L,+2 - \,)Sy (3.24) 
To prove this, start with the formula 

= + (-3i?' - i?2 + V)dl + (-3i2' -E^ + V)'d^ 

- R'" - {VR)' - RR" - R^V, (3.25) 

which implies that {d^ + R){dl + V+){dx - R) = {dx - R){dl + V-){dx + R) is equivalent to 

V± = -R"/R ± 3R' -R^ + C/R. (3.26) 

Now set R := aQ' /Q. Plugging the following identities 



R' = aHl - -^Q'-'). R' = -aP-^QP-\ 
p+l p+l 

p + l 



(3.27) 



into (|3:26]) . we get, for C = 0, 



V± = -a" + -^{a ±{p- l))(a ±{p- 1)/2)QP'\ (3.28) 
p+l 



Hence for a = kj we have 



V± = -kj + -^kj±2kj±i, (3.29) 
■' p + l 

which gives the desired identity (j3.24p . The above proof also shows that Lj-i and are 
the unique choice for the identity to hold with Sj (modulo a constant multiple of Q/Qx, 
which is singular). 

3.5 Variational formulations for eigenvalues of C 

We considered two variational formulations for nonzero eigenvalues of C in general dimen- 
sions in Section 12. 2i Here we present a new variational formulation for 1-D. Define the 
selfadjoint operator 

H := SL+S*. (3.30) 

This is a fourth-order differential operator, with essential spectrum [l,oo). By a direct 
check, we have 

HQ = SL+S*Q = SL+{-2Q^) = 0. 
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Thus Q is an eigenfunction with eigenvalue 0. Since {Q,S*f) = {SQ,f) = for any /, we 
have 

Range S* _L Q. (3.31) 
In particular, since L+|q± is nonnegative for p < 5 by Lemma 12.21 so is H. 

Lemma 3.5 The null space of H is 

N{H) = span{Q,6'P^xQ} , 

where, recall, 5p^ is if p ^ Pc, o,nd 1 if p = Pc- 

Remark. Note that dim N{H) = 1 + 6p^ which is different from dim N {L^^l"^ Lj^L^J^"^) = 
2 + 5p^. We will show below that H and l}'!'^ L^l}!'^ have the same nonzero eigenvalues. 

Proof. If Uf = 0, then L+S*f = -2cQ and S* f = cQi + dQ^ for some c, d G M 
by Lemma 12.11 We have Qi J- Q p = pc = 5. Thus, if p 7^ 5, c = by (j3.3ip . and 
S*{f + fQ) = 0. We conclude / = -f Q. 

When p = 5, we have S*xQ = -Q'^d^iQxQ) = -2Qi. Thus S*{f + fQi + f Q) = 
and / = -fQi - fQ. □ 

Define eigenvalues of H as follows: 

with a suitably normalized minimizer denoted by /j, if it exists. By standard variational 
arguments, if fij < 1, then a minimizer fj exists. By convention, if fi^ is the first of the /Xj's 
to hit 1 (and so fk may not be defined), we set ^ij := 1 for all j > k. 
We can expand Summarv 12.51 to the following. 

Theorem 3.6 (Equivalence) Let n = 1. Let fij be defined as in Summarv \ 2. 5\ and fij he 
defined by (|3.32p . Then fij = jlj. When Hj 7^ and Hj < 1, the eigenf unctions of (j2.19p 
and ()3.32p can be chosen to satisfy 

""i = fj = — SL+Uj. 

Proof. First we establish the equivalence of nonzero eigenvalues. Suppose / = fj is an 
eigenfunction of (|3.32p with eigenvalue /i 7^ 0, then SL^S*f = flf. Let u := S* f ^ and 
apply S* on both sides. By L_ = S*S we get L-L^u = fiu. Thus u is an eigenfunction 
satisfying ()2.14p with fj, = ft. On the other hand, suppose u satisfies L-L+u = with 
/X 7^ 0. Applying SLj^ on both sides and using L_ = S'*S', we get SL^S* SL^u = ^SLj^u, 
i.e., Hf = flf for / = fi'^SL+u. 

Now use Lemmas 12.31 and 13.51 If p G (1, 5), then m = jli = 0, corresponding to Qx and 
Q, and /X2 = > 0. If p = 5, then fii = H2 = f^i = ^-2 = 0, corresponding to Qx, Qi, and 
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Q,xQ, and fi-^ = jl^ = 1. If p G (5,cxd), then = /ii < 0, ^2 = f>'2 = 0, corresponding to 
Qx and Q, and /i3 = ^3 = 1. We have shown fij = fJ-j- □ 

In the following we will make no distinction between and fij. By the minimax 
principle, (13.32P has the following equivalent formulations: 



inf 



{f,Hf) 
sup ,„ 

dim M=j f^M [J,jj 

Here M runs over all linear subspaces of 



sup 



inf 



{f,Hf) 



dimM=j-l /-LA-/ (/, /) 

with the specified dimension. 



(3.33) 



3.6 Estimates of eigenvalues of C 

In this subsection we prove lower and upper bounds for eigenvalues of C, confirming some 
aspects of the numerical computations shown in Figure [TJ Recall that, by Lemma 12.31 the 
first positive fij is fi2 for p £ (l,Pc) and ^3 for p £ \pc,Pmax)- The first theorem concerns 
upper bounds for fii and fi2- 

p 



1 .2 1 



1 .S2 2 



3 ° 

^ - O.J 









M2 












^1 



















^2 



logio(p- 1) 
Figure 6: p vs. fij. 



Theorem 3.7 Suppose n = I and 1 < p < 00. 

(a) If p 7^ 3, then fi2 ^ Cp for some explicitly computable Cp < 1. In particular f2 
exists. 

(b) fii < if and only if p > 5. For any C > 0, we have Hi{p) < —Cp^ for p sufficiently 
large. 

Proof. For part (a), we already know /i2 = for p > 5. Assume p S (1, 5). Considertest 
functions of the form / = SQ^ with A; > 0. / is odd and hence f -L Q, the 0-eigenfunction 
of H. Since H = S'L+S'* and S*S = L_, we have 

^ if,Hf) _ {L_Q\L+L^Q^) 

ifJ) {Q\L^Q^) ■ 

By formulas (|3.5p and (j3.6p . 

L_Q^ = aQ^+P-^ + bQ^, a = -^{k - l){2k + p + I), b=l-k^. 
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Thus 



where 



a = -^—(k + 2p-l)(2k+p-3), d=l-(k+p-l?. 
p + 1 

L+Q'' = cQ'^+P-^ + bQ\ c = -^{k+p){2k -p-l). 

p+l 

(/, Hf) _ a?(Tj^ + a{ad + hc + 6o-) J2 + h{ad + ab + bc)Ji + b^Jo 
if J) aJi + bJo 

Jm= [ Q^''+"'^P-^\x)dx, (m = 0,l,2,3), 



(3.34) 



which are always positive. If A; — > 0"*", then Jm converges to /kQ™^^"^) dx for m > 0, and 
Jg = 0{k^^). The above quotient can be written as 



([33i]) = b' + 



2 . J 



aJi + bJo 
where 

J = a^aJs + a{ad + bc + 6(t) J2 + b{ad + bc)Ji. 
Note that Jra\k=o = ^ /m sech2-(y) dy with sech^"'{y) dy = 2, |, f for m = 

1,2,3, respectively. Also, as A; ^ 0''', a —>■ —1, 6 — > 1, c — > —p, a — > ^'^^^-p^^f-^ , and 
d ^ 1 — [p — 1)^. Direct calculation shows 

hm J = -r—^—T (P + ^)Hp - 3)^ 
k^o+ lb{p — 1) 

Also note 6^ < 1 for k > 0. Thus, if 1 < p < cxd and p ^ 3, then J < and the quotient 
(j3.34p is less than 1 for k sufficiently small. (If p = 3, the sign of J is unclear and (j3.34p 
may not be less than 1.) This proves ;U2 < 1 and provides an upper bound less than 1 for 
fi2- It also implies the existence of f2- This establishes statement (a). 

For statement (b), the fact that ;Ui < if and only if p > 5 is part of Lemma 12.31 We 
now consider the behavior of fii for p large. Fix /c > 1 to be chosen later. As p — > 00, 

= (^)^+"^ • ^ • / (sechx)i^+^™ dx ~ CmP"'-\ 
2 p-l 

with Cm = 2^"*" Jjg( sechx)^™(ix = 2, |, ^ for m = 1, 2, 3, respectively, and 
a~A; — 1, 6 = 1 — k'^, c ~ —p, a ~ 2p, d ~ — p^. 

Thus, by I^M), 

(/, Hf) aaJs + adJ2 1 — k n 

W i^" 
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By choosing k > 1 sufficiently large, we have shown that for any C, ni < —Cp^ for p 
sufficiently large. □ 

The next theorem bounds eigenvalues of C by eigenvalues of and L_ . Recall pj and 
Xj{p) are defined in ()3.4p and (|3.3p . 

Theorem 3.8 (Interlacing of eigenvalues) Fix k > 1 and p G [pk+2, Pk+i) where, re- 
call, pj = Let Xjip) = 1 — + — j{p— 1)]^ he as in (13. 3|) and so \k+i < 1 < 
For the eigenvalues (ij defined by ()3.32p . we have 

A,V(P) < < A|+2(P), {l<j<ky, A^^i(p) < < 1. (3.35) 

In particular, there are K simple eigenvalues ■ ■ ■ , l^K+i in (0, 1) where K = k if fj,k+i < 1 
and K = k — 1 if jik+i = 1- Moreover, K is always 1 when k = 1. Finally, 



JA2A3 (1<P<2), 
A2 (2<^j<5), 



A3A4 (l<p<5/3), 
A3 (5/3 <p<2), 
1 {2<p<oo), 



f^i>-^ip-lf{p-^) (5<p<oo). 
lb 

Remark 3.9 In view of the above lower bounds for /i2 and fi^, we conjecture that 

f^j+i > Xj+iXj+2 {l<p< Pj+2); f^j+i > Xj+i {pj+2 <P< Pj+i)- (3.36) 
This is further confirmed numerically for j = 3,4,5 (see Figure[7]) . Note that limp^p .^-^ 



1 because both Aj+i and fJ-j+i converge to 1. It also seems that ^^^^.'^^^^ has a limit as 

p — > 1+, but it is not clear although we have (|3.35p and Xj = (j — l)(p — 1) + 0{{p — 1)^) 
as p ^ 1+. 

Proof. We first prove the upper bound: For j < k, use the test functions 

5*^2, Sips, S'ljjj+2 

(we cannot use S^i since it is zero). Recall L-tpm = Xmipm- Let a = (02, . . . , aj+2) vary over 
- {0}. By equivalent definition ([333]) . H = SL+S*, L_ = S*S, and the orthogonality 
between the "ipmS, we have 

2 \3 



< sup '" ■ — — - = sup 



m 



< max A2^ = a^+2- 
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Since fJ-j+i < -^^+2 < 1; it is attained at some function, for which the second inequahty 
above cannot be replaced by an equahty sign. Thus /ij+i < '^^4.2- 

For the lower bound of eigenvalues, we use only the special case j = 1 of (13.240 : 

H = SL+S* = SLoS* = S*L3S. (3.37) 
In particular, we have for 1 < p < 3, 

H > S*L2S = S*SS*S = Ll = Ll, (3.38) 

which implies that 

X^j+i < fij+i (1< p < 3) (3.39) 

(and again, equality is impossible). 

For the second eigenvalue fi2, we can get a more precise estimate by using (|3.18p for 
L3 > A3 together with 

Li\q±>X'2, (3.40) 
which follows from spec{Li). Combining these estimates, we have for any f -L Q and p < 5, 

{Hf, f) > A'3(5/, Sf) > X'^X'^if, /), (3.41) 

which implies that fi2 ^ A3A2, i.e., 

[A2 (2<p<5) 
For p > 3, we have is > A2 = — (p — l)(p — 5)/4 and 

L3-L2>-ip-l)ip-3)/2=:-a. (3.43) 
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Hence for any t G [0, 1], we have 

L3 > tL2 - at + (1 - t)X2. (3.44) 

and so for 6 > 0, we have 

{Hf, /) + b{f, /) > {S* {tL2 - at + (1 - t)X2)Sf, f) + /) 

= t\\L^ff - {at - (1 - i)A2)(Li/, /) + b\\ff, ^ ' ^ 

which is nonnegative if 

b>{at-{l-t)X2f/{4t), (3.46) 

whose infimum is attained at t = — A2/(a + A2) = (p — b)/{p — l) G (0, 1) for p > 5. Plugging 
this back in, we obtain the lower bound 

/xi>A2(a + A2) = -^(p-l)^(p-5) (p>5). (3.47) 



We have a similar bound on 113 by using the even-odd decomposition (M) = Lj^ ^ 

-2 



L^^(R). Let ipj, be the eigenfunction of Li and L3 such that 

Litl^j = Xjipj, L-iij = Xjij. (3.48) 

V'j starts from j = 1 and .^j starts with j = 3. They are even for odd j and odd for even j. 
For any even function f ± Q = ipi, Sf is odd and so we have f ± ipi = Q,ijj2 and 5/ _L ^3. 
Hence by spec{L3) and spec(Li), we have 

(Hf, f) = (LsSf, Sf) > X^iSf, Sf) = A4(Li/, /) > A4A3(/, /), (3.49) 

where we denote 

Xr.= h (^<^<^^)' (3.50) 

' \l ipj<p). ^ ' 

Thus the second eigenvalue of H on L^^ is > A4A3. Next for any odd function / ± ■02, we 
have / _L '01, '02, '03- Hence we have 

(Hf, f) = {LsSf, Sf) > X'^iSf, Sf) = A'3(Li/, /) > A'3A4(/, /). (3.51) 

Similarly, every odd function f ± S*^3 satisfies / _L -01 and S'/ ± ^3, ^4, so 

{Hf,f)>XMfJ)- (3.52) 

Hence the second eigenfunction on is > max(A4A3, A5A2) > A4A3. Therefore we have 
> A3A4, i.e., 

'A3A4 (l<p<5/3), 
M3 > < A3 (5/3 <p<2), (3.53) 
.1 (2<P). 

This argument, however, does not yield any useful estimates for the higher fij. □ 
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3.7 Resonance for p = 3 

In the theory of dispersive estimates for the linear Schrodinger evolution, it is important 
to know whether or not the endpoints of the continuous spectrum of the linear operator 
are eigenvalues or resonances. For our £, the endpoints are X = zizi. Resonance here refers 
to a function (p which satisfies the eigenvalue problem locally in space with eigenvalue i or 
—i, but which does not belong to L'^(W^). For dimension n = 1, one requires (j) £ L°°(M). 
(Note for comparison's sake that in one dimension, the operator —(f'/dx^ has a resonance 
- corresponding to the constant function - at the endpoint of its continuous spectrum.) 

Before we made the numerical calculation, we did not expect to see any resonance. 
However, from Figure [H one sees that k = converges to 1 as p — > 3. What does the 
point K = 1 at p = 3 correspond to? A natural conjecture is that it is a resonance or 
an eigenvalue, since the p = 3 case is well-known to be completely integrable and special 
phenomena may occur. 

This is indeed the case since we have the following solution to the eigenvalue prob- 
lem (j2.12p when p = 3, 

\ = i. (3.54) 

It is clear that (f) G L°°(M) but L''(M) for any g < oo. 

Let Up{x) denote the real- valued (and suitably normalized) solution of (12.141) corre- 
sponding to /X = /X2. It is the first component of the eigenfunction of (j2.12p . A natural 
question is: does Up{x) converge in some sense to u^{x) := 1 — Q'^{x) as p ^ 3? Since 
Up — Us is not in L'^{M.) for all q £ [l,oo), it seems natural to measure the convergence in 
the following weighted norm, 

ll/llw := / w(/)2(x) dx, 
JR 

where a weighting operator w is defined by w(/)(x) := f{x) ^^^_^^2 ■ This de-emphasizes the 
value of Up — for x large, and so it should converge to as p goes to 3. This is confirmed 
numerically as follows. 

Let Us := 1 — and 6 := HitsHw In Appendix we will propose a numerical method to 
solve for the eigenpair {A, [up{x), Wp{x)]^} of (j2.12p corresponding to ^2 = — A^. Renor- 
malize Up{x) for p / 3 so that it is real-valued, Up{0) < 0, and ||tip||w = ^- In Figure EJc) 
we plot U3 in a large interval |x| < 130 with 5 = 1.3588. According to the numerical 
method in Appendix, we get ^2.8)^^2.91^3.1 and n3.2 plotted in Figure [8l||a) , (b), (d) and 
(e), respectively. The vertical range is roughly [—1, 1]. In Figure [8l^f)-(j) we plot w{up) for 
p = 2.8,2.9,3,3.1 and 3.2, for |x| < 130 and vertical range [—1,0.5]. 

In Figure [9] we plot p vs. \\up — u^W^ and observe that Up{x) converge to u^^x) in the 
weighted norm ||-||^ as p ^ 3. In the numerical calculation for Figure [9l our increment for 
p is 0.01. 

Remark 3.10 For the operators L^L^ and L^L^, and in general 4-th order operators, it 
seems difficult to exclude the possibility that ^ = 1 is an eigenvalue. Consider the following 
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example. Let H := with p = \/8 — 1- Note —1 is an eigenvalue of when p = \/8 — 1- 

Hence 1 is an eigenvalue of H, at the endpoint of its continuous spectrum. 

It would be interesting to prove the above convergence analytically and characterize the 
leading order behavior near p = 3 as we did in Theorem 12. 61 

4 Excited states with angular momenta 

In this section we consider excited states with angular momenta in M", n > 2. Let k = [n/2], 
the largest integer no larger than n/2. For x = (xi, . . . , Xn) £ M", use polar coordinates rj 
and 9j for each pair X2j-i and X2j, j = 1, . . . ,k. P. L. Lions [20] considers solutions of the 
form 

Q(x) = 0(ri,r2, . ..n,Xn) e^(-i^i+--^-^'=), mj G Z. 

The dependence of (j) in Xn is dropped if n is even. He proves the existence of energy 
minimizing solutions in each such class. 
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For the simplest case n = 2, Q{x) = (j){r) e*"^^ and, by (jl.3p . 4> = '^('^) satisfies 

-</)"- J<A' + ^0 + '^-|<^r-V = O, (r>0). (4.1) 
The natural boundary conditions are 

lim r^™(/)(r) = a, lim r^™^^(/)'(r) = ma, lim (f){r) = 0, (4-2) 

for some q > 0. One can choose (j){r) real-valued. It is shown by laia-Warchall that 
()4.ip - ()4.2p has countably infinite many solutions, denoted by 4'm,k,pif), each has exactly k 
positive zeros. They correspond to "m-equivariant" nonlinear bound states of the form 

Qm,fc,p = 0m,fc,p(Oe""', (A: = 0,1,2,...). (4.3) 

Note that Qm,k,p are radial if and only if m = 0, and the ground state Q = Qofl,p is 
considered in the previous sections. The uniqueness question of (l)m,k,p{f) is not addressed 
in [16j . It is proved for the case A; = in |23j . 

Mizumachi [23j-[26j considered the stability problem for these solutions. He showed 
that 

1. Under m-equivariant perturbations of the form e(r)e™^, Qm,o,p are stable for 1 < p < 
3 and unstable for p > 3; 

2. Under general perturbations, Qm,k,p are unstable for p > 3 for any k; 

3. linear (spectral) instability implies nonlinear instability; (it can be also obtained by 
extending the results in |30J to higher dimensions using the method of pLj); 

4. For fixed p > 1, if m > M{p) is sufficiently large, Qm,k,p are linearly unstable and its 
linearized operator has a positive eigenvalue. 

We are most interested in the last result. Intuitively, for 1 < p < cxd, Qm,k,p should 
be unstable for all (m, k) ^ (0, 0) since they are excited states. Can this be observed 
numerically? It turns out to be true for p away from 1, but false for p near 1. 

In the following, we first describe our numerical methods for k = and next discuss 
their relations. We will discuss our numerical results in the end. We only compute m = 1, 2 
but the same methods work for other m. 

Remark 4.1 Our numerical methods do not apply when A; > 0. Indeed, for m > and 
A; > 0, the radial functions 4'm,k,p{i") are sign-changing and cannot be numerically calculated 
using the method described in the Appendix. In fact, it is an open question whether they 
are unique. Assuming the uniqueness, one needs to develop a new algorithm to compute 
them before one can compute the spectra of jO for Qm,k,p- 
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4.1 Numerical algorithms 

There are two steps in our numerical method: First, compute 4'm,o,pif)- Second, compute 
the spectra of the discretized hnearized operator around Qm,o,p- The second step is more 
involved and we will present three algorithms. 

Step 1. Compute (p{r) = (j)m,o,p{f)- It is energy minimizing among all solutions of 
(|4.ip - (|4.2p for fixed m,p, and it is positive for r > 0. Since our algorithm in the Appendix 
is applicable to all positive (ground state) solutions, we can use it to calculate the discretized 
vector of (p{r) with a small change of the code. 

Step 2. Compute the spectra of the discretized linearized operator. The linearized 
operator C has a slightly different form than (jl.Op because Q = <i>m,o,p{f)e^^^ is no longer 
real. With the same ansatz (ll.6p - (11.7p . the linearized operator C has the form 



Ch = i(^^h-h + ^±^\Q\P~^h + ^|Q|P-3Q2/,j _ (4 4) 

We have developed three algorithms for computing the spectrum of C. 

Algorithm 1. Write Q = cj){r)e'^™'^ = (f)(r) cos{'m6) + i(/)(r) sin(m^). In vector form 
with C acting on [Re/i,Im/i]''', we have 



C 



-A + 1 
A- 1 



+ 



-{p — 1) cos sin — cos —p sin 
p cos^ + sin^ (p — 1) cos sin 



{mO). (4.5) 



h is convenient to use polar coordinates to discretize the operator. We use a two 
dimensional mesh, 

2d mesh: r = : 5^ : Tmax; = Q : 5e : 2tt. (4-6) 

The discretized matrix has size NT by NT with = r^ia^/dr and T = 271/80. We use zero 
boundary condition with Tmax = 

15, 6r = 0.04, and T = 160. 
Although the matrix operator (14. 5p is slightly more complicated than (II. 9p and the mesh 
is 2-dimensional, the same numerical routine can be applied to compute the spectrum of 
the discretized matrix of ()4.5p . The only difference is that the matrix size is much larger. 



Algorithm 2. By restricting the problem to some invariant subspaces of £, as we did 
for the computation of figures HHSl we will reduce the problem to 1-dimension. 

Observe that functions of the form a(r)e*-'^ with a fixed integer j are not preserved by 
C unless j = m, but the following L^-subspaces are invariant under £: 



Xk = AP = |/i(r, 9):h = a{r)e'^'''+^>'' + 6(r)e*^""-'=^''| , < j G N. 

If /c = 0, we drop b{r) and Aq = {/i(r, 6) : h = a(r)e*™'^}. We will compute the spectra of 
C limited to each subspace Xk- Define 
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For A; = 0, with a = oi + ia2 and ai, 02 G M, we have 

£[a(r)e^'"^] = -i [/7oa - Va] e^'"^ = [/7o(a2 - ^i) + ^(02 + iai)]e'"'^ . 
Thus, acting on [01,02]^, C\xo has the matrix form 



Lxo 



Ho + V 
-Ho + V 



For A; > 0, with a = ai + ia2, b = bi + ib2 and ai, 02, ^2 G we have 

= [Hk{a2 - iai) + F(62 + i6i)]e*("+'=)^ + [i?_fc(&2 - ibi) + F(a2 + iaOle*^"-'^)^ 
Thus, acting on [01,02,61,62]^! ^\xk has the matrix form 



To discretize the operator, we use the one-dimensional mesh 











V 


-Hk 





V 








V 







V 





—H-k 






Id mesh: 



0: 6r : r„ 



(4.7) 



The matrix corresponding to Xq has size 2N by 2N. The matrix for Xk with > has 
size 4iV by 4A^. We use zero boundary condition with rmax = 30 and 5r = 0.01. 

Counting multiphcity, the eigenvalues of C is the union of eigenvalues of C\x^, with 
k = 0,1,2,.... 

Algorithm 3. Instead of the form (jl.6p . include the phase e*"*^ in the linearization: 

ip = {(j) + /i)e*'"^^**. Then the linearized operator acting on [Re /i, Im /i]"'' is 



—2m/r'^d0 
-A + 1 + m^/r^ 



-A + 1 + mVr2 - 
-2m/r'^de 



which is invariant on subspaces Zk = {[ai(r),a2(r)]^e*'^^} with integers k. We have 



02 M 



ikd f^ikO T 



ai[r} 
02 (r) 



where 



L 



m,k •" 



-(-Ar + 1 + 



2irnk 



2imk 



9-1 



acting on radial functions. We use the same one-dimensional mesh (14. 7p as in Algorithm 2. 
For every A;, the matrix size is 2N by 2A. We then compute the spectra of Lm,fc foi" each k. 
Counting multiplicity, the eigenvalues of £ is the union of eigenvalues of Lm k with 



k = 0,±1,±2,.... 
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4.2 Properties of these algorithms 

We now address the relation between these algorithms. First note that is essentially 

the sum of and Let us make it more precise and suppose k > 0. The case k = 

is easier. A function h = {ai + m2)(r)e^('"+^)^ + {h + ii62)(r)e^('"-^)^ in Xk C L'^{M.'^) 
can be identified with [01,02,61,62] G = L^jj^(M^; M^). The space X^ is a subspace of 
^radi^^'i *^^) '^^ which we compute the spectrum. The function h can be also identified with 



ai(r) 
02 (r) 



Mr) 
62 (r) 



-ike 



the collection of which form a subspace of © with real components. 

NuUspace of C. The nullspace of C gives a good test of the correctness of our 
numerical results. For k = 0, the 0-eigenfunction iQ of £. corresponds to [0, ^]^e*"*^ in 

p— 1 ^ 

'^''™'^,C), they also provide two 



Xq and [0,(/>]^ in Zq. The generalized eigenfunction Qi = ^^^Q + ^ ' VQ corresponds to 



[Qi,0]Te*™^ in Xo and [^^ + r^',0] ' in Zq. Since Xq C 
(generalized) eigenvectors for Algorithm 1. 
For A; = ±1, the 0-eigenfunctions 



2Qxi = 2(0' cos e-ii; sin 6l)e 



imO 



j(m+l) 



+ (0' + V)e' 



i(m— 1)5 



i(m-l)0 



2Q^2 = 2(<^' sin + cos 6)3'"^^ = i{-(t)' + V')e^("*+^)^ + i((/)' + V)e 
where ■0 = ■m(j)/r, belong to Xi, and correspond to 0-eigenvectors [</>' — ^, 0, + -0, 0]"'' and 



in Zi 



where 



p cos d 
sin ^ 



of Lxi- For Algorithm 3, they correspond to the following vectors 



W+e'^ + W-e 



(p sm U 
ipcosO 



W± 



= -iW+e'^ + iW-e' 



Thus VF+e*^ is a 0-eigenvector of C in Zi, and VF_e 
The generalized eigenfunctions 



is a 0-eigenvector of JC' in Z_i. 



ixiQ = irct) cos Oe'""^ = zr</)e^("^+^)^ + zr^e^^"*-^)^ 



ix2Q = ir (f) sin 6 



r(j)e' 



i(m+i)e 



roe 



,i(m-l)e 



also lie in Xi and correspond to generalized 0-eigenvectors [0, rcf), 0, r(f)]^ and [rcf), 0, — r^, 0]^ 
of Lxi- For Algorithm 3, they correspond to [0,r(j)cos6]^ and [0,r(l)sm9]^ in Zi Z^i. 
By the same consideration as for Q^^ and (5^2 7 their span over C is the same as the span 
of [O,r0]^e*^ G Zi and [O,r0]^e-*^ G Z_i. One can check that 



L 



m,±l 



"0" 






r(f) 


= -2 





(4.8) 
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Thus, the multiplicity of 0-eigenvalue in each of Xq, Zq and Zi is at least 2. The 
multiplicity of 0-eigenvalue on Xi is at least 4. 

Symmetry of spectra. If 



L 



then 



m,—k 



A 



m,—k 



m.k 



A 

-B 



A 



-A 



A " 
-B 



L 



m.k 



A" 

-B 



-A 



-B 



In particular, if A € a{L^^k), then —A G a{Ljn^i;), and A, —A G a{L^^_k). Thus fc) 
itself is symmetric w.r.t. the imaginary axis, and a{Lfn,k) £^nd a{Lfn-k) are symmetric 
w.r.t. the real axis. 

Similarly, one can show that the spectra of Lx^, are symmetric with respect to both real 
and imaginary axes. 

Equivalence of Algorithms 2 and 3. In Algorithm 2, for > 0, we can write 



Lx,= 



HkJ VU ' 
VU H_kJ 



where 



Let 



J 



" 1" 


, u = 


"0 


1" 




"1 0" 





, 1 = 


-1 


1 




1 



M 



I -J 
I J 



I I ' 
J -J 



Noting JU = —UJ, we have 



M-^Lx,M 



aJ + cU -(3 ' 
-P aJ + cU 



1 0" 

10 

10 

1 



a + c 

-a + c 
-(3 

-p 



= P. 








where 



Let 



a 



1 k"^ 



p 



:{Hk - H_k) 



—a + c 

2mk 





P 
a + c 




V. 



L' 



p-^M-^Lx^MP 






P 


a + c 





-p 








a + c 


a + c 








P 





—a + c 


-P 
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In Algorithm 3, L„i^k acts on [A{r), B{r)] . If we write the enlarged matrix of Lm,fc 
acting on [Ke A,lm A,Ke B ,1m B]^ , the matrix is exactly L'. The matrix for L^-k will be 
also L' if it acts on [Re A, — ImA, Rei?, — Imi?]"*^. This amounts to a choice of assigning J 
or — J to the complexification of i. 

More precisely, if L^^^u = An with u = [A, B]^ , then L^k'^u = Xiu. Write A = Ai + iA2 
and B = Bi + iB2 and suppose A; > 0. These two equations are equivalent to 



'a: 




"Re \A 




'-A2 




"-Im \A' 


A2 




\iii\A 


, L' 


Ai 




Re 


Bi 




ReXB 


-B2 




-\m.XB 


B2 




ImXB 




Bi 




ReXB _ 



Adding the second equation multiplied by —i to the first equation, we get 

L'w = Xw, w = [A,-iA,B,-iB]^ . 

Taking conjugation we get L'w = Xw. Thus A and A are eigenvalues of L' , and hence of Lx,,- 
Since Lm,kU = An iff Lm-kU = An, eigenvalues of L^-k also correspond to eigenvalues of 

Counting eigenvalues. Lx^. acts on L^^^(R^,C'^) while Lm,±k act on L^^^(M^, C^). 
The eigenvalues of is the union of eigenvalues of Lm,fc and Lm,-k- For any ball Br on 
the complex plane disjoint from the continuous spectrum Ec = {ir : r S M, |r| > 1}, 

#(c7(LxJ n Br) = #{a{Lm,k) n Bn) + n Br 

which is equal to 2^{a{Lm,k) H Br) if the center of Br is on the real axis. 

Numerical efficiency. Algorithm 1 is 2-dimensional, and thus more expensive to 
compute and less accurate. Both Algorithms 2 and 3 are one-dimensional and more accurate. 

The benefit of Algorithm 3 than Algorithm 2 is that it further decomposes the subspace 
of L^(IR^,C'*) corresponding to Xk to two subspaces. Although its matrix size is only half 
that of Algorithm 2, its components are complex and hence require more storage space. 
Numerically these two algorithms are not very different. 

4.3 Numerical results 

The results of our numerical computations of the spectra of C for m = 1,2 and various 
k and p are shown in Figures fToVllSl As before, we focus on eigenvalues in the square 
{a + bi : \a\ < 1, |6| < 1}. Purely imaginary eigenvalues with modulus greater than 1 corre- 
spond to the continuous spectrum of C, and are discrete due to discretization. 
Let us first describe some simple observations: 

1. The distribution of eigenvalues, see Figures [TOlllH is more complicated and interesting 
than Figures [iHSl There are not only purely imaginary eigenvalues and real eigenvalues 
but also complex eigenvalues, whose existence implies instability. 
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In Figures [T21 - I131 we compare the results obtained from three algorithms. For Algo- 
rithms 2 and 3 the parameter k ranges from to 9. Results from these algorithms 
have high degree of agreement, except when p near 3 and eigenvalues near 0. We will 
discuss this exceptional case in the end. 

For Algorithms 2 and 3, the numerical 0-eigenvalue occurs only when A; = and k = 
±1. It agrees with our discussion in the previous subsection. Their multiplicities also 
match and there is no unaccounted eigenvector. In particular, Ng{C) has dimension 
6 if p 7^ 3 and 8 if p = 3, the same as the case of ground states. We also numerically 
verified the nullspace, for example, the discrete version of (|4.8|) is correct. 

As p increases, two pairs of purely imaginary eigenvalues may collide away from 0, 
and then split to a quadruple of complex eigenvalues which are neither real nor purely 
imaginary. For m = 1, this bifurcation phenomenon appears three times before p = 
1.55 and there are 3 complex quadruples for p > 1.55. For m = 2, it occurs five times 
before p = 1.5 and there are 5 quadruples for p > 1.5. These complex eigenvalues 
seem to move away from the imaginary axis as p increases further. 

As p increases to 3, (by Algorithms 2 and 3), a pair of purely imaginary eigenvalues 
from 0-th subspace collide at and then split to a pair of real eigenvalues as p increases 
further. This is the same picture as in the ground state case in section 3. Indeed, 
Mizumachi [23] proves that Qm,o,p are stable in the 0-th subspace if p < 3 and unstable 
if p > 3. Thus p = 3 is a bifurcation point. Also note that when p = 3 the NLS (jl.ip 
has conformal invariance and explicit blow-up solutions can be found as in the ground 
state case. 

In Figures [1^151 we observe the bifurcation more closely. For m = 1, the bifurcation 
occurs when {k,p,X) equal 

(0, 3, 0), (1, 1.52765 -0.436i), (2, 1.0165, -0.016i), (3, 1.3495, -0.219i). 

For m = 2, the bifurcation occurs when {k,p,X) equal 

(0, 3, 0), (1, 1.357, -0.180i), (2, 1.007, -0.027i), (3, 1.0245, -0.035i) 

and 

(4, 1.0455, -0.045i), (5, 1.3955, -0M7i). 

Due to the existence of complex eigenvalues for m = 1, 2 and p > 1.02, Qm,o,p is spec- 
trally unstable for these parameters. However, all these complex eigenvalues bifurcate 
from some discrete eigenvalues ±bi with |6| < 1 and p > 1.008. Our computation for 
both m = 1, 2 and 

p = £- 0.001, ^ = 1,2,3,... ,8, (up to 15 if m = 1) 
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does not find any complex eigenvalues. This suggests that the two excited states 
(j)ifl,p{r)e^^ and 4'2,o,p{r)e'^'^^ are linearly stable when p is sufficiently close to 1. It is 
possible that the numerical error increases enormously as p — > due to the artificial 
boundary condition, since the spectrum is approaching to the continuous one for p = 1. 
This has to be verified analytically in the future. 

We finally discuss the exceptional case when p is near 3 for eigenvalues near 0. In this 
case Algorithm 1 produces an quadruple of complex eigenvalues ±0.0849 it 0.0836i, and the 
0-eigenvalue has multiplicity 4. We expect to see larger errors from Algorithm 1 but the 
error in this exceptional case is much larger. It is related to the large size of a Jordon block 
for the eigenvalue. As discussed in the previous subsection, the nullspace is at least 6 
dimensional. The analysis in section [231 suggests that (we do not claim a proof), as p goes 
to the bifurcation exponent pc = 3 from below, a pair of imaginary eigenvalues merges into 
the Jordon block containing the eigenfunctions iQ and Qi, and the Jordon block becomes 
size 4. As is well-known in matrix analysis (see [9l p. 324], if a matrix contains a 

Jordon block of size i, the computed eigenvalues corresponding to that block have errors of 
order e^/^, where e is the sum of the machine zero, the truncation error from discretization, 
and the perturbation (from varying p). Since 6r = 0.04 for Algorithm 1 and the truncation 
error of a central difference scheme for has order 0(5^), the error for the zero eigenvalue 
near p = 3 could be 

(,52)1/4 ^ 

In contrast, for other bifurcation points on the imaginary axis, the Jordon block at the 
bifurcation exponent is of size 2 and the error is of order (5^)^/^ = 0.04. In practice, the 
error is smaller due to cancellation and the numerical results by Algorithm 1 do not differ 
too much from those by Algorithms 2 and 3. Also note that numerically the eigenspace 
has dimension 4, accounting for Q^^ and ixjQ. The complex quadruple correspond to iQ, 
Qi and the joining pair of nonzero eigenvalues. 

Appendix: Numerical method 

In this section we describe a numerical method to compute the spectrum of the linear 
operator C defined by (jl.Sh for p > 1 and space dimension n > 1. There are two main steps 
in this method. First, we will solve the nonlinear problem (II. 3p for Q: we will discretize 
it into a nonlinear algebraic equation, and then solve it by an iterative method. Second, 
we will compute the spectrum of C: we will discretize the operator £ into a large-scale 
linear algebraic eigenvalue problem and then use implicitly restarted Arnoldi methods to 
deal with this problem. 

Hereafter, we use the bold face letters or symbols to denote a matrix or a vector. For 
A € R^^^^, q = {qi, . . . , qn)"^ S M.^ , q® = q o • • • o q denotes the p-time Hadamard 
product of q, and [[q| := diag(q) the diagonal matrix of q. 

Step I. We first discretize equation ()1.3p into a nonlinear algebraic equation and 
consider it on an n-dimensional ball 17 = {x G : |x| < R,R ^ M}. We rewrite the Laplace 
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operator —A in the polar coordinate system with a Dirichlet boundary condition. Based 
on the recently proposed discretization scheme |19j . the standard central finite difference 
method, we discretize — Aq(x) into 

Aq = A[gi,...,gjv]^, AgM^x^, (A.l) 

where q is an approximation of the function Q{x). The matrix A is irreducible and 
diagonally-dominant with positive diagonal entries. The discretization of the nonlinear 
equation (jl.3p can now be formulated as the following nonlinear algebraic equation, 

Aq + q-q® =0. (A. 2) 

We introduce an iterative algorithm [15] to solve ()A.2p : 

Aq,+i+q,+i =q®, (A.3) 

where q^+i and q^ are the unknown and known discrete values of the function Q{x), re- 
spectively. The iterative algorithm is shown below. 



Iterative Algorithm for Solving <5(x). 

Step Let j = 0. 

Choose an initial solution Hn > and let qn = ii~'°ii . 

Step 1 Solve the equation (|A.3p . then obtain qj+i- 

Step 2 Let a^+i = and normalize qj+i to obtain q^+i = aj+iq^+i. 

Step 3 If (convergent) then 

1 

Output the scaled solution (Qj+i)p-i qj_|_i. Stop. 

else 

Let j := j + L 
Goto Step 1. 



If the components of qo are nonnegative, this property is preserved by each iteration qj, 
and hence also by the limit vector if it exists (see [151 Theorem 3.1]). The convergence of 
a subsequence of this iteration method to a nonzero vector is proved in [15[ Theorem 2.1]. 
Although the convergence of the entire sequence is not proved, it is observed numerically to 
be very robust. See Chen-Zhou-Ni [5j for a survey on numerically solving nonlinear elliptic 
equations. 



Step II. Next we discretize the operator £ of (jl.lOp into a linear algebraic eigenvalue 
problem: 

(A.4) 



u 


= A 


u 






w 




w 
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where 



L = 





A-I+bq®l 



A + I-[q®] 




7 = p — 1, u = (ui, . . . , Un)~^ G w = {wi, . . . ,W]\f)~^ G M^, and q is the output of 
the previous step, and satisfies the equation in ()A.2p . We use ARPACK [2T] in MATLAB 
version 6.5 to deal with the linear algebraic eigenvalue problem ()A.4p and obtain eigenvalues 
A of L near the origin for p > 1 and space dimension n > 1. Furthermore, the eigenvectors 
of L can be also produced. 

The Step II above can in principle be used to compute all eigenfunctions in L^(M"). 
However, in producing Figures [IHSl we look for eigenfunctions of the form 4>{r)e^'"^^ . These 
problems can be reformulated as 1-D eigenvalue problems for (/>(r), which can be computed 
using the same algorithm and MATLAB code. This dimensional reduction saves a lot of 
computation time and memory. Even with this dimensional reduction, and applying an 
algorithm for sparse matrices, the computation is still very heavy, and we cannot compute 
all eigenvalues in one step. We can only compute a portion of them each time. 
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Figure 10: Spectra of C in for m = 1 and < k < 24 as p 
1.06, 1.1, 1.4, 1.55, 1.8, 2.3, 3, 3.1, 3.5 computed by Algorithm 3. 
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Figure 11: Spectra of C in for m = 2 and < k < 28 as p 
1.06, 1.1, 1.3, 1.4, 1.8, 2.2, 2.6, 3, 3.5 computed by Algorithm 3. 
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Figure 12: Spectra of C in for m = 1 and various p = 1.6,2.1,3,3.2. Point "•" denotes 
the spectra computed by Algorithm 1 and the others symbols denote the spectra computed 
by Algorithms 2 and 3. 
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Figure 13: Spectra of C in for m = 2 and various p = 1.6,2.1,3,3.2. Point "•" denotes 
the spectra computed by Algorithm 1 and the others symbols denote the spectra computed 
by Algorithms 2 and 3. 
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Figure 14: Bifurcation diagrams of Lm,k for m = 1 and < /c < 3. 
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Figure 15: Bifurcation diagrams of Lm,k for m = 2 and < A; < 5. 



